class particula
{
   private:

   public:
       ///Posicion
	   double pos[3];
	   ///Posicion anterior
	   double posOld[3];
	   ///Velocidad
	   double vel[3];
	   ///Fuerza
	   double force[3];
	   double drg[3];
	   double drv[3];
	   ///Posicion de la particula en determinado momento
	   double rstart[3];
	   ///Celda en la que se encuentra la particula
	   int celda;
	   ///Coordenadas de la celda
	   int celdaxyz[3];
       ///Numero de vecinos
	   int nvecinos;
       ///Lista con indice de vecinos
	   int lstVecinos[1000];
	   ///Radio
	   double radio;
	   double mass;

	   particula()
	   {
            radio=1.0;
	   }

	   particula(double rr[3])
	   {
		   for (int i=0;i<3;i++)
		   {
			   pos[i]=rr[i];
			   posOld[i]=rr[i];
			   rstart[i]=rr[i];
		   }
		   radio=1.0;

	   }

	   ///posicion de la particula
	   void setPosition(double rr[3]);
	   ///Velocidad de la particula
	   void setVelocities(double vv[3]);
	   ///Fuerza de la particula
	   void setForces(double ff[3]);
};

	   void particula::setPosition(double rr[3])
	   {
		   for (int i=0;i<3;i++)
		   {
			   pos[i]=rr[i];
			   posOld[i]=rr[i];
			   rstart[i]=rr[i];
		   }
	   }

	   void particula::setVelocities(double vv[3])
	   {
		   for (int i=0;i<3;i++)
			   vel[i]=vv[i];
	   }

	   void particula::setForces(double ff[3])
	   {
		   for(int i=0;i<3;i++)
		   {
			   force[i]=ff[i];
		   }
	   }


